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NONLINEAR  ANALYSIS  OF  A  MAGNICON  OUTPUT  CAVITY 


I.  Introduction 

The  generation  of  the  microwaves  required  to  power  future  high-gradient  radio- 
frequency  (rf )  accelerators  will  require  the  use  of  efficient  and  powerful  sources  that 
operate  at  higher  frequencies  than  conventional  S-band  klystrons.  Possible  rf  sources 
to  meet  this  need  include  the  relativistic  klystron,  the  gyroklystron  and  the  magni- 
con.  In  addition,  there  are  other  potential  applications  for  such  sources,  including 
microwave  power  beaming.  This  paper  presents  an  analysis  of  the  electron-wave 
interaction  in  the  output  cavity  of  a  magnicon.  The  magnicon1-6  is  an  advanced 
version  of  the  gyrocon6  and  employs  a  scanning  beam  that  is  obtained  by  the  pas¬ 
sage  of  a  magnetized  pencil  beam  from  the  electron  gun  through  a  deflection  system. 
The  deflection  system  consists  of  an  input  cavity  and  one  or  more  passive  cavities, 
separated  by  drift  tubes,  with  the  entire  system  immersed  in  an  axial  magnetic 
field,  B0.  The  cavities  support  a  rotating  TM110  mode  with  a  frequency  that  is 
~  1/2  the  gyrofrequency,  wc  =  |e|i?o/7w*c.  Here,  e  is  the  charge  and  m  is  the  mass 
of  an  electron,  7  is  the  relativistic  factor,  and  c  is  the  vacuum  speed  of  light.  The 
purpose  of  the  deflection  system  is  to  spin  the  beam  to  high  transverse  momentum; 
i.e.,  a  =  t >x/vx  >  1.  Here,  Uj.  and  vt  are  the  velocity  components  transverse  to  and 
along  the  z  axis.  After  passing  through  the  deflection  system,  the  beam  transverse 
momentum  is  used  to  drive  a  gyrotron-like  interaction  in  the  output  cavity.  The 
entry  point  of  the  electrons  in  the  output  cavity  rotates  in  space  about  the  cavity 
axis  at  the  drive  frequency.  In  the  frequency-doubling  version,  the  output  cavity 
supports  a  rotating  TM210  mode  with  frequency  u>  «  u>e.  Each  electron  arrives  in 
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the  output  cavity  ideally  gyrophased  for  optimum  energy  transfer  to  the  rf  field. 
Since  the  electrons  entering  the  output  cavity  are  almost  completely  phase-bunched 
and  rotate  in  synchronism  with  the  TMjio  wave,  the  transverse  efficiency  may  be 
extremely  high. 

At  the  Naval  Research  Laboratory  (NRL)  a  program  to  develop  a  high-gain, 
second-harmonic  magnicon  in  the  X-band  regime  is  presently  under  way.3  The  design 
calls  for  the  generation  of  50  MW  at  11.4  GHz  and  50%  efficiency,  using  a  200  A, 
1/2  MV  electron  beam  produced  by  a  cold-cathode  diode  on  the  NRL  Long-Pulse 
Accelerator  Facility.7 

This  paper  presents  an  investigation  of  the  output  cavity  of  a  magnicon  and  is  a 
follow-up  of  our  earlier  analysis  of  the  deflection  system.6  The  equations  of  motion 
for  the  electrons  and  the  wave  equation  for  the  rf  field  are  derived.  The  derivation 
closely  follows  the  analysis  of  the  quasi-optical  gyrotron  presented  in  Refs.  8-10.  In 
particular,  it  is  assumed  that  the  build-up  of  the  rf  field  in  the  output  cavity  takes 
place  on  a  time  scale  that  is  long  compared  to  the  transit  time  of  electrons  through 
the  cavity.  Scaled  variables  are  defined  and  the  equations  are  analyzed  in  detail 
to  examine  the  efficiency  of  the  device.  Results  from  time-dependent  and  steady- 
state  simulations  of  the  output  cavity  are  presented.  Point  design  parameters  for  a 
magnicon  operating  in  X-band  are  given.  A  major  simplification  in  this  analysis  is 
the  use  of  the  ideal  TM2io  mode  structure  in  the  output  cavity,  neglecting  the  effects 
of  field  leakage  and  fringing  fields  in  the  vicinity  of  the  beam  tunnel.  Additionally, 
the  effects  of  space  charge  Eire  neglected  in  the  analysis. 
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II.  Formulation 


The  vector  potential  of  the  electromagnetic  field  in  the  cavity  is  assumed  to  be  a 
superposition  of  discrete  Fourier  harmonics 

A(r,t)  =  ^2  ~ txp(- iput)  +  c.c.,  (1) 

p*o  1 

where  A^p\  the  amplitude  of  the  pth  Fourier  component,  may  be  decomposed  in 
terms  of  the  cavity  modes  as  follows 


^(P)  =  5Z  ^w^WPnmr/aJcos^/iJexp^n^),  (2) 

n,mji 


where  r,  <f>,  z  denote  the  cylindrical  coordinates,  is  a  complex- valued  amplitude 

which  is  a  slowly-varying  function  of  time,  t,  Jn  is  the  ordinary  Bessel  function  of 
the  first  kind  of  order  n,  is  the  mth  zero  of  J„,  a  is  the  radius  and  L  is  the 
length  of  the  cavity  and  c.c.  stands  for  complex  conjugate.  Figure  1  shows  a  cross- 
sectional  view  of  the  electric  and  magnetic  field  lines.  The  whole  structure  rotates 
azimuthally  at  the  drive  frequency,  w/2.  The  vector  potential  in  Eq.  (1)  is  a  solution 
of  the  wave  equation  in  the  form 


V2AX 


1  d2Ax  _  4; rl  ft 
c*  &t2  c  Loodio 

x  J  d2rxo  J  d3^  /(rioj  Po, *o)^*(rJ.o>  Po> *o)^[r  —  r(rXo >  Po?  t,  to)].  (3) 


where  /  is  the  beam  current  and  a  tilde  over  a  variable  indicates  that  the  variable 
is  a  function  of  t  and  the  transverse  coordinates  rxo  and  momenta  p0  at  the  entry 
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time  to  into  the  cavity.  The  right  hand  side  of  Eq.  (3)  is  obtained  by  writing  the 
electron  distribution  function  in  the  form 

/(r,p,t)  -  !  dt0 

|e|  J—oo 

x  J  d3r±o  J d3p0/(rjLo,Po,<o)£[r  -  r(rxo,p0,Mo)]%  -  p(rj.o,Po,Mo)], 

where  /(rxo,  Poj*o)  is  the  distribution  function  at  the  entrance  to  the  cavity.  The 
distribution  function  /(r,  p,  t)  is,  by  construction,  a  solution  of  the  Vlasov  equation. 


Herein,  the  interest  is  in  the  excitation  of  the  TM2i0  mode,  represented  by  the 
Fourier  amplitude  j4jjq,  with  no  z  dependence  within  the  cavity.  Making  use  of  the 
Fourier- Bessel  expansion  of  5[t  —  r(t,<o)|  a-nd  substituting  Eqs.  (1)  and  (2)  into  Eq. 
(3),  one  obtains 


]Texp(— ipwf) 


2ipu  d  ,  (pW  i4«o  _  4//c  fL  dz 

.  c*  dt  +  \  c*  VJ  *  [aJ1(p7i)]i  Jo  L 


X  f  dtnVj6(z 

J — OO 


*)Jz{PnTla)  exp(-2»<£), 


(4) 


where  the  cutoff  wavenumber  is  defined  by  ke  =  pz\/a  —  wfc.  In  writing  Eq. 
(4),  the  right-hand  side  has  been  simplified  by  assuming  that  there  is  no  spread 
in  the  transverse  coordinates  or  momenta  of  the  electrons  as  they  enter  the  cavity. 
The  cavity  radius  a  is  determined  by  the  boundary  condition  J2(P2i)  —  0.  Since 
P21  =  5.136,  the  radius  of  the  output  cavity  is  2.145  cm  for  a  frequency  of  11.42 
GHz. 

Assuming  that  the  Fourier  amplitude  j42iq  varies  slowly  on  the  time  scale  2tt/uj, 
and  performing  the  integration  (2i r/u»)_1  fo*'”  dt  exp(iut)  on  Eq.  (4),  the  wave 
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equation  for  Aji],  takes  the  form 

2iw  d  (u1  2\1  cEo(i)  _  4/ fc  rL  dz 
c3  dt  »  c2  c  1 1  tu>  [<x«/i(p2i  )]*  Jo  L 

Xlo  2tt j <j)  Jl  {P3i^[T(<o,2),io]/a)}«cp{-2t^[r(to,«),<o]  +  (5) 

where  E0(t)  =  wA$0/c  is  the  slowly  varying  complex- valued  electric  field  amplitude 
and  r(to,z)  =  to  +  /o  dz’/vz(t o,z')  is  the  time  measured  along  the  electron  orbit.  As 
discussed  in  Ref.  11,  in  going  from  Eq.  (4)  to  Eq.  (5),  the  order  of  integration  has 
been  interchanged  according  to  the  prescription  /0J^“  dt  dt0  — *  So"^  dt0  /!„,  dt. 

Next,  it  is  necessary  to  write  down  the  equations  of  motion  of  the  electrons.  In 
terms  of  the  independent  variable  2,  these  equations  take  the  form 


dvx 

dz 

\e\£Zv  (\t\B*vy  |e|fltf\ 

7T7ic2  "m  \  7 me  vz  7mc  J  ’ 

(6) 

dVy 

dz 

lel£«r  f\e\B* 

7 me2  9  \  7 me  7 me  vz  J  ’ 

(7) 

dvz 

\e\£z  /c2  \  (\e\Byvm 

(8) 

dz 

7 rod  \d,  *)  \~fmcvz  ymcvz)' 

where  £,  =  —(l/c)dAM/dt  is  the  rf  electric  field,  (Bm,Bv)  =  Vx(i,e,)  is  the  rf 
magnetic  field  and  e,  is  a  unit  vector  along  the  z  axis.  Since  the  electrons  interact 
with  fields  that  are  independent  of  the  z  coordinate,  it  follows  that  the  canonical 
momentum  along  the  z  axis  is  conserved;  i.e., 


Pz  =  7  mvz  — 


|e|  A, 


—  const. 


(9) 


5 


To  proceed,  it  is  convenient  to  transform  to  guiding  center  variables  in  the  guide 
field  Bo,  defined  by 


Vg  —  1>xCOS$, 

Vy  =  Esin'S, 

*  —  X  +  psin$, 
y  =  Y  —  pcos 


where  $  is  the  gyroangle,  p  =  vx/ujc  is  the  gyroradius,  u>c  =  |e|B0/7mc  is  the 
relativistic  gyrofrequency,  and  (X,Y)  is  the  guiding-center  coordinate  in  the  (x,y)- 
plane.  The  equations  of  motion  of  the  guiding-center  are 


dX  B. 
dz  ~  B0’ 

dY_  ^  B± 
dz  Bo 


(10) 

(11) 


The  equations  of  motion  for  the  other  variables  may  be  simplified  by  approximating 
the  Bessel  function  as  follows: 


•MP2i*7«)exp(-2*^)  » 

x  [BJ  exp(-2*0)  -  p2exp(-2t$)  -f  2*pBexp[-»($  +  0)]j  , 


where  Bexp(*0)  =  X  +  iY  defines  the  guiding  center  radius,  R,  and  the  guiding 

center  angle,  0.  The  equations  of  motion  of  the  electrons  along  the  cavity  may  now 

be  written  in  the  form 

dux  E 

7/7o 


R2  sin^  -  2pBsin(d'  —  6)  +  p2  sm(i>  —  26) 


6 


+ 


€ 


(12) 


[iZsin(y»  —  6)  —  psin(V’  —  20)j, 

7/7o 

— - Lfi3  sin  V>  -  2p-Rsin(V>  -  fl)  -f  p3  sin(V>  —  20)1 

o(  7/7o  P*o««  1  J 

—  — ~~ — [iisin(^  —  9)  -  /5sin(V>  —  20)],  (13) 

7/7o 

~  =  -^-~~-L-—[Rcob^-0)-pcos(^-2$)),  (14) 

7/7o  ux  7/7o«i 

where  £  =  z/L  is  the  distance  along  the  cavity  scaled  to  the  cavity  length,  R  — 
PnR/2a  and  p  —  p2ip/2a  are  the  scaled  guiding- center  radius  and  gyroradius,  re¬ 
spectively,  i*x  —  vi./v*o  and  uM  =  ux/v*o  are  the  components  of  the  electron  velocity 
scaled  to  the  z  component  of  the  velocity  at  the  entrance  to  the  output  cavity,  u*oi 
Pm  —  vzo/c>  =  \e\BoLfa0mcvM  is  the  scaled  strength  of  the  guide  magnetic  field, 
7o  is  the  relativistic  factor  at  the  entrance  to  the  cavity,  E0{t)  =  *|J50(*)|  exp[tjS(t)] 
defines  the  phase  0{i)  of  the  rf  field,  e  =  \eE0\L/'y0mcvt0  is  the  scaled  strength  of 
the  rf  electric  field  and 

t  =  (vL/vM)f‘dOu.(C)-0-2O  o.  (15) 

Figure  2  indicates  the  definition  of  the  single  9  =  $  —  0  +  ?r/2  in  terms  of  the 
guiding  center  angle  0  and  gyroangle  4>.  These  angles  rotate  at  the  drive  frequency 
(w/2)  such  that  at  the  injection  point  to  the  cavity  9  does  not  vary  with  time;  i.e., 
0(C,to)  =  0O  +  w<0/2,  $(C,*o)  =  HO  +  wf0/2,  where  ©0  is  a  constant  and  $(C) 
describes  the  variation  of  the  gyroangle  along  the  cavity. 

When  integrating  the  equations  of  motion,  Eqs.  (10)-(15),  along  the  cavity,  the 
rf  field  is  taken  to  be  a  constant.  This  is  a  consequence  of  the  assumption  that  the 


7 


temporal  ^volution  of  the  rf  amplitude,  Eo(t),  takes  place  on  a  time-scale  that  is 
mu«*h  longer  than  the  electron  transit  time  (ss  L/v, )  through  the  cavity.8"10  Hence, 
the  use  of  partial  derivatives  in  Eqs.  (12)-(14). 

To  complete  the  set  of  governing  equations,  it  is  necessary  to  recast  the  wave 
equation,  Eq.  (5),  in  scaled  form.  We  are  interested  in  a  system  that  is  phase 
locked  by  a  prebunched  electron  beam.10  It  is  assumed  that  the  prebunched  beam 
induces  an  rf  field  in  the  cavity  at  a  frequency  w  =  ckc  +  Aw,  where  Aw  is  the 
frequency  shift  with  respect  to  the  cold  cavity  frequency.  Intrcducing  the  cavity 
quality  factor,  Q,  into  Eq.  (5),  and  making  use  of  the  small  argument  expansion  of 
the  Bessel  function,  the  wave  equation  reduces  to 

+  =  I  £  d([R3sitni>-2pRsin(T/>-0)  +  p2sin{tl>-20)]  ,  (16) 

e  ^Au»  —  =  ij  d(  (f2J cos^  —  2pRcos(it>  —  0)  +  P1  cos(V>  —  20)]  ,  (17) 

where  /  =  IQL/  is  the  scaled  beam  current,  IA  =  (mc3/M)7o/?*o 

is  the  Alfven  current,  r  =  ckct/Q  is  the  scaled  time  and  Aw  =  Q / ckc  is  the 
scaled  frequency  mismatch.  It  is  important  to  note  that  the  integration  over  the 
entry  time  <0  in  Eq.  (5)  can  be  trivially  performed  and  does  not  appear  in  Eqs. 
(16)  and  (17).  This  is  due  to  the  special  nature  of  the  magnicon,  in  which  electrons 
enter  the  output  cavity  prebunched,  and  also  due  to  our  assumption  that  there 
is  no  spread  in  the  initial  transverse  coordinates  or  momenta.  Thus,  to  study  the 
complete  evolution  of  the  system,  it  is  only  necessary  to  6tudy  the  passage  of  a  single 
electron  through  the  cavity.  This  is  the  basis  of  the  simulation  results  presented  in 
Secs.  Ill  and  IV.  We  shall  address  the  issue  of  non-ideal  beam  effects  in  Secs.  V . 
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The  efficiency  is  defined  as  rj  =  (70  —  7/)/(lo  —  1),  where  7/  is  the  relativistic 
factor  at  the  cavity  exit.  Making  use  of  Eqs.  (12),  (13)  and  (16),  it  follows  that,  in 
steady  state, 


&o7o  ** 

V  2(7o-1)/' 

It  should  be  remarked  that  the  expression  in  Eq.  (18)  is  simply  a  restatement  of 
the  balance  between  the  power  extracted  from  the  electron  beam  and  the  power 
dissipated  in  the  finite- Q  output  cavity. 


In  studying  the  efficiency  of  conversion  of  electron  energy  into  rf  energy  as  an 
electron  traverses  the  output  cavity,  it  is  useful  to  define  the  detuning,  A,  that 
measures  the  deviation  of  the  gyrofrequency  from  the  rf  frequency: 


_  (Wc)I  , 

A  =  'AT'i°- 


(19) 


In  common  with  the  gyrotron,  for  zero  detuning  the  efficiency  should  be  small  and 
significant  conversion  will  be  achieved  only  for  positive  detuning.13  An  estimate  for 
the  optimum  detuning  may  be  obtained  by  choosing  the  magnitude  of  the  guide 
field  so  as  to  maintain  gyroresonance  on  the  average,  assuming  that  the  electrons 
are  brought  to  a  complete  rest  at  the  end  of  the  interaction.  That  is,  B '%*  — 
(mcw/je |)(7o  +  l)/2.  Making  use  of  Eq.  (19),  the  optimum  detuning,  Aopt,  is  given 
by 


Aopt  _  {v/c)L  70  ~  1 

PzO  27o 


(20) 
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III.  Single-Electron  Steady- State  Simulation;  Re¬ 
duced  Equations 

The  purpose  of  this  section  is  to  present  steady-state  simulation  results  obtained 
from  numerical  solution  of  the  reduced  equations  derived  in  the  previous  section. 
Neglecting  guiding- center  drifts,  Eqs.  (12)-(17),  with  8 /dr  set  equal  to  zero,  consti¬ 
tute  a  simplified  set  of  equations  that  can  be  solved  to  rapidly  scan  the  parameter 
space.  Since  these  equations  are  expressed  in  seeded  form,  the  simulation  results 
may  be  used  to  obtain  point  designs  for  a  variety  of  physical  systems. 

Equations  (10)  and  (11)  represent  the  drift  of  the  electron  guiding-center.  The 
right-hard  sides  of  Eqs.  (10)  and  (11)  are  on  the  order  of  10  %  or  less  for  the 
examples  described  in  detail  here.  This,  along  with  the  close  agreement  between 
steady-state  simulation  results  based  on  the  reduced  equations,  and  those  based  on 
the  full  Lorentz  equations  in  Sec.  IV,  justifies  the  neglect  of  guiding-center  drifts. 

All  the  simulation  results  in  this  paper  assume  that  on  entering  the  output 
cavity,  the  guiding  center  of  an  electron  is  located  off  the  cavity  axis  at  a  distance 
equal  to  the  gyroradius.  This  assumption  is  roughly  borne  out  by  our  previous 
investigation  of  the  electron  beam  as  it  is  spun  up  in  the  deflection  cavities.® 

The  input  parameters  for  the  simulation  are  the  beam  energy  at  injection,  the 
scaled  cavity  length,  fceL,  and  the  detuning,  A.  Equation  (20)  is  used  to  obtain 
an  order-of-magnitude  estimate  for  the  detuning  that  leads  to  high  efficiency.  The 
angle  6,  which  is  defined  following  Eq.  (15)  and  in  Fig.  2,  is  related  to  the  gyroangle 
and  the  guiding-center  angle  at  the  entrance  to  the  cavity.  This  angle  is  also  an 
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input  parameter,  being  determined  in  part  by  the  length  of  the  drift  tube  leading 
to  the  output  cavity.  For  a  given  value  of  the  scaled  electric  field,  e,  the  scaled 
current,  I,  and  the  efficiency  are  determined  from  Eqs.  (16)  and  (18),  respectively. 
The  plots  presented  in  this  section  are  obtained  by  scanning  over  all  rf  phase  values, 
— ir  <  /?  <  7r,  to  determine  the  maximum  efficiency  and  the  corresponding  values  of 
the  scaled  frequency  shift,  the  scaled  current,  and  the  rf  phase.  To  leading  order, 
the  approximation  u >  ~  ckc  is  employed  in  Eq.  (15)  and  the  frequency  shift  is  then 
simply  obtained  from  Eq.  (17).  The  last  input  parameter  is  the  initial  beam  a. 
Two  examples  will  be  discussed  in  detail,  corresponding  to  the  initial  values  a  =  1 
and  a  =  1.5. 

Case  (i)  :  a  =  1 

Figure  3  shows  the  results  of  the  single-particle,  steady-state  simulation  of  the 
reduced  equations  for  a  =  1.  The  parameters  for  this  simulation  are  listed  in 
Table  l,  the  three  different  detunings  corresponding  to  the  three  curves  in  each 
plot.  Figures  3(a),  (b),  (c),  and  (d)  show  the  efficiency,  17,  the  scaled  current,  I,  the 
scaled  frequency  shift,  Aw,  and  V’o,  respectively.  The  latter  is  the  value  of  the  phase 
angle  1/)  at  the  entrance  to  the  cavity.  From  Eq.  (15),  it  follows  that  V’o  =  — /?  — 20o, 
thus  relating  Fig.  3(d)  to  the  rf  phase,  /?.  The  abscissa  in  each  plot  is  the  scaled 
electric  field. 

On  examining  Fig.  3,  it  is  interesting  to  note  that  significant  conversion  effi¬ 
ciency  (>  50%)  may  be  obtained  with  a  beam  a  as  small  as  1.  It  is  also  interesting 
to  note  the  discontinuity  in  some  of  the  curves  at  high  values  of  scaled  electric  field. 
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This  discontinuity  results  from  a  new  local  efficiency  maximum  (as  a  function  of 
phase)  becoming  the  global  maxi  mum.  (Our  plots  of  maximum  efficiency  as  a  func¬ 
tion  of  electric  field  do  not  reveal  the  variation  of  the  efficiency  with  rf  phase;  see 
Ref.  10.) 

It  should  be  noted  that  the  existence  of  steady-state  solutions  does  not  guarantee 
their  accessibility  or  their  stability  in  an  experiment  that  must  grow  from  rf  noise. 
In  order  to  investigate  the  accessibility  of  the  high-efficiency  solutions  that  the 
scanning  code  locates,  we  will  make  use  of  a  time-dependent  code,  as  discussed  in 
Sec.  IV. 

Case  (ii)  :  a  =  1.5 

Figure  4  shows  the  results  of  the  single-particle,  steady-state  simulation  of  the 
reduced  equations  for  a  —  1.5.  The  parameters  for  this  simulation  are  listed  in 
Table  II,  the  three  different  detunings  corresponding  to  the  three  curves  in  each 
plot.  Figures  4(a),  (b),  (c)  and  (d)  show  the  efficiency,  tj,  the  scaled  current,  /,  the 
scaled  frequency  shift,  Aw,  and  V’oj  respectively.  The  general  behavior  of  the  plots 
is  similar  to  that  in  Fig.  3,  although  no  discontinuity  is  apparent  in  this  case.  Note 
that,  as  expected,  the  efficiency  can  be  quite  high,  reaching  up  to  70  %  for  the  case 
with  A  =  5.788. 
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IV.  Single- Electron  Time-Dependent  Simulation; 
Maxwell-Lorentz  Equations 

In  Sec.  Ill  we  presented  a  panoramic  view  of  the  magnicon  output  cavity  based 
on  the  highly  simplified  equations  derived  in  Sec.  II.  Figures  3  and  4  permit 
one  to  design  an  output  cavity  by  translating  the  scaled  variables  according  to  the 
constraints  of  any  particular  experiment. 

It  is  pointed  out  in  Sec.  Ill  that  the  maximum  efficiency  steady-state  pre¬ 
dicted  by  a  time-independent  code  may  not  be  accessible  or  stable  in  an  experiment 
in  which  the  fields  build  up  from  noise.  To  address  this  issue  one  can  employ  a 
time- dependent  code  to  determine  which  one  of  the  optimum-efficiency  final  states 
indicated  in  Figs.  3  and  4  are  true  steady  states  of  the  system.  The  time-dependent 
simulation  results  presented  here  are  obtained  from  the  full  set  of  Maxwell-Lorentz 
equations,  Eqs.  (5)-(8),  following  the  motion  of  a  single  electron  through  the  cavity 
and  studying  the  build-up  of  the  rf  field  over  a  much  longer  time  scale. 

The  parameters  for  the  time-dependent  simulation  are  chosen  by  locating  a 
high  efficiency  point  on  the  solid  curve  in  Fig.  3  or  Fig.  4,  and  reading  off  the 
corresponding  value  for  the  scaled  current,  and  the  frequency  shift.  An  important 
criterion  in  choosing  an  operating  point  is  to  minimize  the  electric  field,  consistent 
with  an  acceptable  efficiency,  in  order  to  avoid  breakdown  problems.  In  the  exam¬ 
ples  considered  in  this  paper,  the  maximum  tolerable  electric  field  is  taken  to  be 
ss  300  kV /cm.  A  further  constraint  on  the  parameters  relates  to  the  length  of  the 
cavity.  Conservation  of  canonical  momentum,  Eq.  (9),  implies  that  the  kinetic  z- 
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momentum,  px  =  “fmvty  can  vary  along  the  cavity.  We  find  that  there  is  an  optimal 
cavity  length  for  which  the  efficiency  is  a  maximum  and  p,  is  a  minimum  as  the 
electron  exits  the  cavity.  This  point  will  be  clarified  presently.  The  current  /  to  be 
used  in  the  time-dependent  code  iB  determined  by  taking  a  reasonable  value  for  the 
cavity  quality  factor,  Q  (<  1000).  The  choice  for  the  cavity  length  and  the  detuning 
determine  the  value  of  the  guide  field  through  Eq.  (19).  Having  chosen  the  param¬ 
eters,  the  time-dependent  code  is  used  to  search  for  the  maximum  efficiency  by 
fine-tuning  the  frequency  shift  around  the  value  determined  from  Fig.  3(c)  or  Fig. 
4(c).  It  must  be  emphasized  that  the  time-dependent  simulation  results  are  based 
on  the  full,  unapproximated  set  of  Maxwell-Lorentz  equations  and  may  therefore 
differ  to  some  extent  from  those  based  on  the  reduced  equations  of  motion. 

Case  (i):  a  =  1 

Table  III  lists  the  parameters  for  the  time-dependent  simulation  of  the  case  with 
initial  beam  a  =  1,  A  =  3.776  and  /  =  6.5.  Figures  5(a),  (b),  and  (c)  show  plots 
of  the  electric  field  amplitude,  the  rf  phase,  and  the  efficiency  as  functions  of  time. 
Figure  5  shows  that  the  system  settles  into  a  steady  Btate  after  a  transient  that  lasts 
w  60  ns.  Figure  5(a)  shows  that  the  electric  field  builds  up  to  about  305  kV/cm, 
Fig.  5(b)  shows  that  the  phase  settles  to  an  asymptotic  value  1.1  rad,  and  Fig.  5(c) 
shows  the  final  efficiency  in  this  field  is  about  52  %. 

Figures  6  (a),  (b),  and  (c)  show  plots  of  the  axial  momentum  (normalized  to 
its  initial  value),  efficiency,  and  a  for  an  electron  traversing  the  output  cavity  in 
the  final  steady-state  field.  The  modulation  of  the  electron  momentum,  pz,  is  seen 
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to  be  correlated  with  the  modulation  on  the  efficiency  curve.  The  scale-length  for 
the  modulation  is  «  2 ir(3t/ke.  The  length  of  the  cavity  has  been  chosen  to  be  4.75 
cm  so  that,  on  exiting  the  cavity,  pt  is  a  minimum,  and  therefore  the  efficiency 
is  a  maximum.  In  this  sense,  it  appears  that  there  is  an  optimal  length  for  the 
output  cavity,  which  may  be  fine  tuned  to  obtain  the  maximum  efficiency.  Figure 
6(c)  shows  that  the  value  of  a  declines  as  the  electron  spins  down  on  traversing  the 
cavity. 

In  terms  of  the  scaled  variables,  the  final  state  of  the  time-dependent  code  is 
compared  with  the  corresponding  point  on  Fig.  3  in  Table  IV.  This  table  shows 
close  agreement  between  the  results  from  the  steady-state  simulation  (labeled  SSS) 
of  reduced  equations  and  the  time-dependent  simulation  (labeled  TDS)  of  the  full 
Maxwell-Lorentz  equations. 

It  must  be  pointed  out  that  the  steady-state  simulations  are  deficient  in  a  very 
important  respect.  Namely,  the  questions  of  accessibility  and  stability  of  the  pre¬ 
dicted  final  states  of  the  cavity,  shown  in  Fig.  3,  remain  unresolved.  As  a  case  in 
point,  we  compare  in  Table  V  the  results  of  the  time-dependent  and  steady-state 
codes  at  the  point  at  which  the  scaled  current  /  =  17.13  and  the  detuning  A  =  4.776, 
which  lies  at  the  extreme  right  of  Fig.  3  on  the  dotted  curve,  and  ostensibly  corre¬ 
sponds  to  the  highest  efficiency.  Examination  of  Table  V,  however,  reveals  that  for 
this  point  the  time-dependent  code  asymptotes  to  a  final  state  with  a  much  smaller 
efficiency.  This  example  highlights  the  fact  that  not  all  of  the  optimum-efficiency 
states  in  Fig.  3,  or  Fig.  4,  are  accessible  and  stable  final  states  of  the  rf  field  that 
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is  built  up  from  noise  in  the  output  cavity. 

Case  (ii):  a  =  1.5 

Table  VI  lists  the  parameters  for  the  time-dependent  simulation  of  the  case  with 
initial  beam  a  =  1.5,  A  =  4.788  and  I  =  6.75.  Figures  7(a),  (b)  and  (c)  show  plots 
of  the  electric  field  amplitude,  the  rf  phase  and  the  efficiency  as  functions  of  time. 
This  figure  shows  that  the  system  settles  into  a  steady  state  after  a  transient  that 
lasts  «  25  ns.  Figure  7(a)  shows  that  the  electric  field  builds  up  to  about  230 
kV/cm,  Fig.  7(b)  shows  that  the  phase  settles  to  an  asymptotic  value  -1.2  rad,  and 
Fig.  7(c)  shows  the  final  efficiency  in  this  field  is  about  63  %. 

Figures  8(a),  (b)  and  (c)  show  plots  of  the  axial  momentum  (normalized  to  its 
initial  value),  efficiency,  and  a  for  an  electron  traversing  the  output  cavity  in  the 
final  steady-state  field.  The  length  of  the  cavity  has  been  chosen  to  be  6.2  cm  so  that 
p ,  is  a  minimum  at  the  end  of  the  cavity  and  therefore  the  efficiency  is  a  maximum. 
Figure  8(c)  shows  that  the  electron  spins  down  on  traversing  cavity,  although  not 
to  the  same  extent  as  that  in  Case  (i). 

In  terms  of  the  scaled  variables,  the  final  state  of  the  time-dependent  code  is 
compared  with  the  corresponding  point  on  Fig.  4  in  Table  VII.  Examination  of  Table 
VII  shows  close  agreement  between  the  results  from  the  steady-state  simulation  of 
reduced  equations  and  the  time-dependent  simulation  of  the  full  Maxwell- Lorentz 
equations.  A  comparison  has  also  been  made  between  the  simulations  for  I  =  10.9 
on  the  solid  curve  in  Fig.  4;  close  agreement  between  the  results  is  obtained. 
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V.  Multi-Electron  Steady-State  Simulation;  Lorentz 
Equations 

On  the  basis  of  the  single-electron  simulation  of  the  preceding  sections,  point  designs 
for  operation  of  a  magnicon  output  cavity,  with  a  =  1  or  a  =  1.5,  have  been  ob¬ 
tained.  The  purpose  of  this  section  is  to  assess  the  role  of  beam  quality,  as  measured 
by  spreads  in  the  input  parameters  of  the  beam,  on  the  efficiency.  The  numerical 
results  presented  in  this  section  are  obtained  by  solving  the  Lorentz  equations  of 
motion  for  a  collection  of  electrons  that  represent  spreads  in  the  initial  values  of  a, 
the  guiding  center  radius,  R,  or  7.  In  all  cases,  the  electron  distribution  is  taken  to 
be  a  “top-hat”,  with  initially  uniform  loading  of  the  electrons.  This  suffices  to  give 
an  order-of-magnitude  estimate  of  the  effect  of  beam  quality  on  the  efficiency. 

The  parameters  for  the  multi-electron,  steady-state  simulation  are  listed  in  Ta¬ 
bles  III  and  VI.  The  asymptotic  rf  phase  obtained  from  the  time-dependent  simula¬ 
tion  is  used  as  an  input  parameter  to  the  steady-state  code.  The  electron  equations 
of  motion  are  integrated  along  the  cavity  for  several  field  amplitudes  to  pin-point 
the  rf  steady-state  that  is  in  power  balance  with  the  given  electron  current.  In  the 
figures,  the  spreads  are  given  in  terms  of  the  root-mean-square  (rms)  deviation  from 
the  mean  of  the  corresponding  quantity.  The  abscissa  in  the  figures  represent  the 
ratio  of  the  rms  to  the  mean  value,  expressed  as  a  percentage. 

Before  proceeding,  it  is  important  to  remark  on  a  basic  premise  of  the  steady- 
state  simulation  results  of  this  section.  As  the  spread  in  a  beam  parameter  increases 
from  zero,  the  rf  field  at  which  power  balance  is  obtained  varies.  For  a  sufficiently 
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large  value  of  the  spread,  however,  there  is,  in  some  cases,  a  sudden  jump  in  the 
rf  field  amplitude  required  for  power  balance.  In  the  figures  to  be  described,  the 
maximum  value  of  the  spread  is  limited  to  lie  below  the  point  at  which  the  sudden 
jump  takes  place.  Nevertheless,  it  cannot  be  guaranteed  that  a  multi-electron, 
time- dependent  code  would  lead  to  results  identical  to  those  presented  here. 

Case  (i):  a  =  1 

Figures  9  (a),  (b),  and  (c)  show  the  decline  in  the  efficiency  due  to  increasing 
spreads  in  a,  guiding  center  radius,  R,  and  7,  respectively.  The  rapid  decrease  of 
efficiency  with  energy  spread  is  specially  noteworthy.  Obviously  this  implies  that 
the  design  of  the  deflection  system  must  be  particularly  cognizant  of  this  sensitive 
dependence  on  energy  spread. 

Case  (ii):  a  =  1.5 

Figures  10  (a),  (b),  and  (c)  show  the  decline  in  the  efficiency  due  to  increasing 
spreads  in  a,  guiding  center  radius,  12,  and  7,  respectively.  Once  again,  one  finds 
that  the  optimum  efficiency  is  extremely  sensitive  to  a  spread  in  7.  Additionally, 
comparison  of  Figs.  9  and  10  shows  that  the  beam  with  the  higher  a  is  relatively 
mor»*  sensitive  to  spreads  in  the  beam  parameters. 
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VI.  Summary  and  Conclusions 


The  electron- wave  interaction  in  the  output  cavity  of  a  frequency-doubling,  magnetic- 
field-immersed  magnicon  has  been  analyzed  in  a  series  of  simulation  studies.  In  this 
presentation  the  output  cavity  is  considered  as  an  entity  that  is  separate  from  the 
deflection  system,  with  the  electron  beam  parameters  corresponding  to  what  is 
expected  on  the  basis  of  our  earlier  studies  of  the  deflection  system.8 

Single-electron,  steady-state  simulation  of  reduced  equations  of  motion  were 
used  to  rapidly  scan  the  parameter  space  and  locate  a  desirable  operating  point.  A 
single-electron,  time-dependent  code  was  then  employed  to  ensure  that  the  chosen 
final  state  was  an  accessible  and  stable  operating  point  of  an  amplifier  wherein  the 
rf  grew  from  noise.  Finally,  a  multi-electron,  steady-state  code  was  used  to  study 
the  sensitivity  of  the  point  design  to  spreads  in  the  beam  parameters. 

In  conclusion,  it  is  useful  to  recall  the  salient  points  of  this  work.  First,  as 
indicated  in  Figs.  6  and  8,  for  an  ideal  beam,  it  is  possible  to  choose  the  cavity 
length  so  as  to  convert  not  only  the  transverse  momentum  but  also  part  of  the  axial 
momentum  into  rf  field  energy.  Second,  based  on  the  runs  made  with  the  time- 
dependent  code,  we  have  found  that  some  of  the  final  states  determined  by  the 
steady-state  code  are  accessible  and  stable,  and  some  are  not.  Third,  an  efficient 
(>  50%)  final  state  with  a  =  1  is  achievable.  Finally,  the  efficiency  of  the  output 
cavity  is  most  sensitive  t,o  energy  spread  on  the  beam. 
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Table  captions 


Table  I:  Parameters  for  steady-state  simulation  of  a  magnicon  amplifier  output  cavity 
for  three  detunings.  Initial  beam  a  =  1. 

Table  II:  Parameters  for  steady-state  simulation  of  a  magnicon  amplifier  output  cavity 
for  three  detunings.  Initial  beam  a  =  1.5. 

Table  III:  Parameters  for  time-dependent  simulation  of  an  X-band  magnicon  amplifier 
output  cavity.  Initial  beam  a  =  1. 

Table  IV:  Comparison  of  single-particle  time-dependent  simulation  (TDS)  and  single¬ 
particle  steady-state  simulation  (SSS)  of  reduced  equations.  Beam  a  =  1, 
detuning  A  =  3.776  and  scaled  current  I  =  6.5. 

Table  V:  Comparison  of  single-particle  time-dependent  simulation  (TDS)  and  single¬ 
particle  steady-state  simulation  (SSS)  of  reduced  equations.  Beam  a  =  1, 
detuning  A  =  4.776  and  scaled  current  I  —  17.13. 

Table  VI:  Parameters  for  time-dependent  simulation  of  an  X-band  magnicon  amplifier 
output  cavity.  Initial  beam  a  =  1.5 

Table  VII:  Comparison  of  single-particle  time-dependent  simulation  (TDS)  and  single¬ 
particle  steady-state  simulation  (SSS)  of  reduced  equations.  Beam  a  —  1.5, 
detuning  A  =  4.788  and  scaled  current  /  =  6.75. 


Table  I 


Voltage 

KL 

11.37  cm 

Beam  a 

DHBHjli 

Detuning  At 

2.776 

Detuning  A  3 

3.776 

Detuning  A3 

4.776 

e 

rc/2 
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Table  II 


Voltage 

500  kV 

keL 

14.85  cm 

Beam  a 

Detuning  Aj 

3.788 

Detuning  Aj 

4.788 

Detuning  A3 

5.788 

6 

DHHHI 

I 

: 
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Table  III 


Frequency 

11.43 

GHz 

Voltage 

500 

kV 

Current  I 

0901! 

Cavity  Radius 

2.145 

cm 

Cavity  Length 

4.75 

cm 

Cavity  Quality  Factor  Q 

200 

Beam  a 

1 

Magnetic  Field  Bo 

6.455 

kG 

Detuning  A 

3.776  | 

Gyroangle  $ 

0  j 

Guiding  Center  Radius  R 

3.2 

mm 

Frequency  Shift  Au>/2rr 

6 

MHz 
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Table  IV 


Ea 

SSS 

e 

2.34 

Au> 

0.11 

■Hi 

-0 

-1.1 

-0.9 

■ 

51.9  % 

48.4  % 
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Table  V 


TDS 

sss 

6 

0.7 

4 

Aw 

0.475 

0.475 

-0 

2.9 

-1.7 

EH 

1.8  % 

57.6  % 
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Table  VI 


Frequency 

11.4138 

GHz 

Voltage 

HEBI 

Current  I 

Cavity  Radius 

2.145 

cm 

Cavity  Length 

6.2 

cm 

Cavity  Quality  Factor  Q 

130 

Beam  o 

1.5 

Magnetic  Field  B0 

6.85 

kG 

Detuning  A 

4.788 

Gyroangle  § 

7f/2 

Guiding  Center  Radius  R 

3.5 

mm 

Frequency  Shift  Aw/27r 

-10.2 

MHz 
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Table  VII 


€ 

2.95 

2.9 

AC) 

-0.12 

-0.12 

-0 

EES 

naanai 

tm 

63  % 

60.6  % 
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X 


Figure  2:  Projection  of  gyro-orbit  onto  transverse  plane.  Electron  is  at  (®,y)  and  guid¬ 
ing  center  is  at  (X,  Y).  (i2, 0)  denotes  guiding  center  radius  and  angle,  and 
(p,  $)  denotes  gyroradius  and  gyroangle.  Angle  0  =  $  -  0  +  tt/2  is  fixed  for 
electrons  entering  output  cavity. 
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Figure  3:  Results  from  single-electron,  steady-state  simulation  of  reduced  equations  for 

initial  beam  a  =  1  and  three  detunings,  Ai  =  2.776  ( - ),  A2  =  3.776  ( - ), 

A3  =  4.776  ( . ).  (a)  Efficiency,  77;  (b)  Scaled  current,  /;  (c)  Scaled  frequency 

shift,  Aw;  (d)  ip0.  Abscissa  is  scaled  electric  field,  c. 


Figure  4:  Results  from  single-electron,  steady-state  simulation  of  reduced  equations  for 

initial  beam  a  =  1.5  and  three  detunings,  Aj  =  3.788  ( - ),  A2  =  4.788 

( - ),  A3  =  5.788  ( . ).  (a)  Efficiency,  tj ;  (b)  Scaled  current,  /;  (c)  Scaled 

frequency  shift,  Aw;  (d)  V'o-  Abscissa  is  seeded  electric  field,  e. 
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Pigure  5:  Results  from  single-electron,  time-dependent  simulation  of  X-band  magnicon 
with  initial  beam  a  =  1  and  detuning  A  —  3.776.  (a)  Electric  field  amplitude, 
|£o|:  (b)  RF  phase,  —  /?;  (c)  Efficiency,  tj.  Abscissa  is  time. 
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Z  (cm) 

Figure  6:  Results  from  single-electron,  time-dependent  simulation  of  X-band  magnicon 

with  initial  beam  a  —  1  and  detuning  A  =  3.776.  (a)  Axial  kinetic  momentum 
normalized  to  initial  value,  pM/pt0]  (b)  Efficiency,  tj;  (c)  a.  Abscissa  is  distance 

along  cavity.  These  plots  correspond  to  motion  of  an  electron  through  cavity 


in  final  steady  state. 
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Figure  7:  Results  from  single-electron,  time- dependent  simulation  of  X-band  magnicon 
with  initial  beam  a  =  1.5  and  detuning  A  =  4.788.  (a)  Electric  field  ampli¬ 
tude,  |220|;  (b)  RF  phase,  -/?;  (c)  Efficiency,  17.  Abscissa  is  time. 
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Figure  8:  Results  from  single  electron,  time-dependent  simulation  of  X-band  magnicon 
with  initial  beam  a  =  1.5  and  detuning  A  =  4.788.  (a)  Axial  kinetic  mo¬ 
mentum  normalized  to  initial  value,  (b)  Efficiency,  rj;  (c)  a.  Abscissa 

is  distance  along  cavity.  These  plots  correspond  to  motion  of  am  electron 
through  cavity  in  final  steady  state. 
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Figure  10:  Results  from  multi-electron,  steady-state  simulation  with  initial  beam  a  =  1.5. 

Efficiency  versus  relative  rms  deviations  :  (a)  <  (a-  <  a  >)J  >,/J  /  <  a  >, 
(b)  <  (R-  <  R  >)*  >l/i  /<*>,( c)  <  (7-  <  7  >)’  >,/S  /  <  7  >•  Here; 

<>  indicates  an  average  over  electron  distribution. 


